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ABSTRACT 

The  radio-controlled  model  (RCM)  has  been  used 
extensively  for  model  tests  since  1970.  An  assumption  that 
has  been  made  is  that  the  angles  measured  by  the  RCM 
gyros  are  the  Euler  angles  for  the  Euler  coordinate 
transformation.  A  recent  study  of  the  gyros  indicates  that 
this  assumption  is  incorrect  not  only  for  the  RCM,  but  for 
most  model  testing  instrumentation  packages  used  at  DTRC. 
(The  Euler  angle  assumption  can  be  a  good  assumption  as 
long  as  the  combination  of  pitch  and  roll  angles  are  small.) 
This  report  presents  a  set  of  transformation  equations  that 
can  be  used  to  convert  RCM  data  to  the  desired  Euler 
angles. 


ADMINISTRATIVE  INFORMATION 

This  effort  was  funded  under  the  Stability,  Control  and  Maneuvering 
Characteristics  program  of  the  SSN  751  Submarine  tests  performed  at  the  David 
Taylor  Research  Center  (DTRC).  This  effort  was  performed  at  DTRC  under 
Work  Unit  1564-147. 


INTRODUCTION 

There  are  a  variety  of  free-running  model  instrumentation  packages  at 
DTRC  that  are  used  to  make  test  measurements  on  models  during  maneuvers. 
These  measurements  are  then  used  to  predict  the  corresponding  full-scale 
responses.  As  part  of  the  instrumentation,  gyros  are  generally  used  to  measure 
the  yaw,  pitch  and  roll  attitude  angles  of  the  model.  These  attitude  angles  are 
then  used  as  predicted  full-scale  attitude  angles  or  as  input  data  for  further 
analysis  using  the  equations  of  motion.  A  recent  study  of  the  radio-controlled 
model  (RCM)  gyros  revealed  that  the  measurements  made  by  these  gyros  are  not 
the  Euler  attitude  angles  as  presumed.  To  measure  the  Euler  angles,  a  stable- 
table  gyro  system,  similar  to  the  gyros  used  on  full-scale  vehicles,  would  be 
required.  Due  to  size  and  money  constraints,  all  of  the  model  experiments 
conducted  at  DTRC  use  a  vertical  gyro  to  measure  pitch  and  roll,  and  a  free 
gyro  to  measure  yaw.  A  vertical  gyro  uses  an  erection  system  to  maintain  the 
spin  axis  in  a  vertical  position;  a  free  gyro  uses  a  caging  mechanism  to  set  the 
initial  direction  of  the  spin  axis.  The  desired  Euler  angles  can  be  computed  by 
processing  the  instrumentation  measurements  with  the  transformation  equations 
derived  in  this  report. 

Vector  analysis,  in  combination  with  coordinate  transformations,  was  used 
to  derive  equations  to  transform  the  RCM  measurements  to  Euler  angles. 
Appendix  A  gives  a  detailed  derivation  of  the  transformation  equations. 

DISCUSSION 

All  displacement  gyro  measurement  systems  rely  on  the  fact  that,  in  the 
absence  of  external  forces  and  moments,  a  rotating  mass  will  maintain  the  same 
axis  of  rotation  (spin  axis)  relative  to  a  fixed  coordinate  system.  An  example  of 
a  vertical  gyro  is  given  in  Fig.  1.  The  gyro  is  constructed  with  two  sets  of 
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Fig.  1.  Gyro  construction  showing  rotating  mass,  gimbals,  and 
angle  measurement  devices. 
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gimbals  that  isolate  the  rotating  mass  from  external  forces  and  moments  caused 
by  changes  in  the  attitude  of  the  gyro  case.  The  displacement  angles  are 
measured  by  using  angular  measurement  devices  mounted  on  the  gimbals  and  on 
the  case  of  the  gyro. 

Three  constructs  are  necessary  to  define  the  angle  that  a  gyro  will  measure, 
that  is,  the  angle  of 

•  the  reference  plane, 

•  the  measurement  plane,  and 

•  the  position  vector. 

Since  the  spin  axis  of  the  rotating  mass  remains  fixed  in  space,  the  spin  axis  can 
be  thought  of  as  a  vector  that  is  perpendicular  to  a  plane  in  the  inertial  (or  fixed) 
coordinate  system.  This  plane  is  the  reference  plane.  The  measurement  device 
(usually  a  potentiometer  or  synchro)  is  designed  to  measure  angular  deflections 
in  the  plane  of  the  device  (or  body).  This  plane  is  the  measurement  plane.  The 
last  construct  necessary  to  define  the  measured  angle  is  the  position  vector, 
which  extends  from  the  center  of  the  measurement  device  to  the  null  position. 
The  measured  angle  is  the  angle  between  the  position  vector  and  the  intersection 
of  the  reference  plane  and  the  measurement  plance.  For  example: 

A  vertical  gyro  that  is  being  used  to  measure  the  model  roll  on  the  outer 
gimbal  potentiometer  has  the  following  constructs: 

•  The  spin  axis  of  the  rotating  mass  is  oriented  in  the  vertical  direction, 
implying  that  the  reference  plane  is  the  horizontal  plane. 

•  The  outer  gimbal  potentiometer  is  attached  directly  to  the  case  which 
means  that  the  measurement  plane  is  the  body’s  Y-Z  plane. 

•  The  position  vector  is  a  vector  from  the  center  of  the  measurement 
device  to  the  null  position  of  this  device. 

When  the  gyro  case  is  displaced  in  roll  (rotation  around  the  outer-gimbal  axis), 
the  outer  gimbal  potentiometer  arm  will  be  positioned  at  the  intersection  of  the 
body’s  Y-Z  plane  and  the  horizontal  (X-Z)  plane.  The  angle  measured  will  be 
the  angular  displacement  between  the  null  position  and  this  new  position.  The 
attitude  angle  definitions  used  in  the  David  Taylor  Research  Center  (DTRC) 
submarine  equations  of  motions  are  given  in  Fig.  2. 

The  purpose  of  the  gyro  measurements  is  to  provide  information  to  relate 
vectors  in  one  coordinate  (fixed  or  body)  system  (such  as  gravity  acting  on  the 
mass  of  the  vehicle)  to  corresponding  vectors  in  the  other  coordinate  (body  or 
fixed)  system.  The  fixed  and  body  vector  quantities  (e.g.,  forces  acting  on  the 
center  of  gravity)  are  related  by  means  of  a  coordinate  transformation  matrix. 
The  Euler  transformation  matrix  for  transforming  fixed  coordinate  system 
measurements  to  the  body  coordinate  system  is  obtained  by  rotating  the  fixed 
coordinate  system  into  the  body  coordinate  system  using  three  angular  rotations. 
The  first  rotation  (see  Fig.  3)  involves  a  rotation  around  the  fixed  Z-axis  until 
the  X-axis  of  the  intermediate  system  is  in  a  position  that  is  directly  under  the 
X-axis  of  the  body  system.  This  angle,  is  the  yaw  angle  of  the  Euler  system. 
The  next  rotation  is  around  the  Y-axis  of  the  intermediate  coordinate  system 
created  by  the  above  rotation  (see  Fig.  3).  During  this  second  rotation,  the  X- 
axis  of  the  intermediate  system  rotates  into  the  X-axis  of  the  body’s  coordinate 
system.  The  second  rotation  angle,  0,  is  the  pitch  of  the  Euler  system.  The  last 
rotation  involves  a  rotation  around  the  body’s  X-axis  (see  Fig.  3).  During  the 
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Fig.  3.  Rotations  used  to  derive  the  transformation  from  the  fixed 
coordinate  system  to  the  model  coordinate  system. 
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final  rotation,  the  second  intermediate  Z-axis  rotates  into  the  body  Z-axis.  This 
angle,  is  the  roll  angle  used  in  the  Euler  system. 

This  sequence  of  rotations  transforms  fixed  (or  inertial)  vectors  into  body 
(or  model)  frame  vectors.  The  inverse  of  this  transformation  (i.e.,  the  Euler  to 
fixed  coordinate  system  transformation)  is  defined  in  Eq.  1.* 

cosif»cos0  ( -  sinipcos^  +  costf»sini^sin6)  (sim^sintp  +  cosi^cosipsinS)  ig 
sim^cosfl  (cost#^cos4  +  sin<j>sintf»sin0)  (sim^sin6cos<f>  -  sin'^cosifi)  jg  (1) 
-  sin0  cos6sin<|i  cosScos^  kg 

The  ip,  jp,  and  kp  are  the  fixed  coordinate  system  unit  vectors  and  the  ig,  jg, 
and  kg  are  the  body  coordinate  system  unit  vectors.  The  variables  ifj,  6,  iji  are 
the  Euler  angles  for  yaw,  pitch,  and  roll,  respectively. 

Note  that  the  coordinate  transformations  being  discussed  consist  of  a 
sequence  of  rotations.  The  individual  matrices  that  describe  each  of  these 
rotations  are  unitary  matrices;  the  inverse  is  equal  to  the  transpose.  The 
transformation  matrix  is  obtained  by  multiplying  the  individual  rotation  matrices 
together.  The  properties  of  the  matrix  transpositions  imply  that  the  inverse  of  a 
transformation  matrix  is  also  a  unitary  matrix;  the  transpose  is  the  inverse.  See 
Appendix  A  for  more  details  on  this  convenient  property. 

Equation  (1)  gives  rise  to  the  following  relationships  between  the  time 
derivatives  of  the  Euler  angles  and  the  components  of  the  body  centered  angular 
rates  measured  using  a  rate  gyro  package.  The  p,  q,  and  r,  given  below,  are  the 
components  of  the  body  angular  velocities  resolved  along  the  body’s  X-,  Y-,  and 
Z-axes,  respectively.  The  variables  4/,  0,  and  ^  are  time  derivatives  of  the  Euler 
angles.  Equations  (2)  through  (4)  express  the  body  angular  rates  in  terms  of  the 
Euler  angular  rates. 

p  =  ^-i^»sin0, 

« 

q  =  0cos'^  +  4»cos0sin'^,  and 

r  =  -0sin^  +  4)cos0cos4. 

Equations  (5)  through  (7)  express  the  Euler  angular  rates  in  terms  of  the  body 
angular  rates.  These  equations  were  obtained  by  solving  Eqs.  (2)  through  (4)  for 
0,  and  u'- 

4  =  p  +  tan0(r  cos<^  +  q  sin^) 

• 

0  =  q  cosi^  -  r  sin^ 

4;  =  sec0(r  cos<ti  +  q  sin"^) 

The  above  sequence  of  rotations  can  be  expressed  in  terms  of  reference 
planes,  measurement  planes,  and  position  vectors  as  shown  in  Table  1. 

•Goldstein,  H.,  Classical  Mechanics,  Addison-Wesiey  Publishing  Co.  (1965).  The  transformation 
was  derived  by  applying  the  methods  presented  in  this  reference  to  the  DTRC  submarine 
coordinate  system  (see  Fig.  2). 


(5) 

(6) 
(7) 


(2) 

(3) 

(4) 
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Table  1.  Position  vectors,  measurement  planes,  and  reference 
planes  for  Euler  angle  measurements. 


Measurement 

Position  Vector 

Measurement 

Plane 

Reference 

Plane 

Yaw 

Vertical  projection 
of  body’s  X-axis 
on  fixed  X-Y  plane 

Fixed  X-Y  plane 

Fixed  X-Z  plane 

Pitch 

Body’s  X-axis 

Vertical  plane 
passing  through 
body’s  X-axis 

Fixed  X-Y  plane 

Roll 

Body’s  Y-axis 

Body’s  Z-Y  plane 

Fixed  X-Y  plane 

The  instrumentation  difficulty  involved  in  measuring  the  Euler  angles  is 
related  to  the  measurement  plane  definition.  To  understand  this  difficulty,  a 
discussion  of  the  measurement  planes  for  displacement  gyros  is  needed. 

The  outer  gimbal  measurement  device,  as  mentioned  earlier,  is  attached  to 
the  gyro  case.  This  makes  the  measurement  plane  a  plane  in  the  body’s 
coordinate  system,  which  implies  that  the  angle  measured  by  this  potentiometer 
is  a  body-centered  measurement.  The  inner  gimbal  measurement  device  is 
mounted  in  such  a  fashion  that  the  spin  axis  will  always  lie  in  the  measurement 
plane.  In  the  case  of  a  vertical  gyro,  this  measurement  plane  will  always  be 
vertical.  As  discussed,  the  Euler  pitch  and  roll  can  be  obtained  directly  from  a 
vertical  gyro  with  the  inner  gimbal  measuring  pitch  and  the  outer  gimbal 
measuring  roll.  The  measurement  plane  for  yaw  is  the  fixed  X-Y  plane.  This 
requires  that  the  readout  device  be  maintained  in  a  horizontal  attitude.  The  only 
way  to  maintain  this  measurement  plane  would  be  to  use  a  gyro  mounted  on  a 
stable  table.  As  mentioned  earlier,  cost,  size  and  weight  considerations  make  this 
measurement  impossible  to  obtain  in  a  direct  fashion  for  model  tests.  The 
normal  instrumentation  used  to  measure  yaw  is  a  free  gyro  using  the  outer 
gimbal  for  the  readout  device.  This  outer  gimbal  yaw  measurement  can  be 
transformed  to  the  desired  Euler  yaw  by  using  Eq.  (8). 

^J  =  tan  “  ’  [(tanifi2cos0  -i-  sin<|>sin0)/cosij>]  (8) 

Where  ip,  6,  and  ^  are  Euler  yaw,  pitch,  and  roll,  respectively,  and  ip2  is 
measured  yaw  angle. 

This  is  the  normal  instrumentation  complement  used  at  DTRC.  The  error 
associated  with  ignoring  the  correction  of  Eq.  8  can  be  seen  from  two  examples. 
If  the  measured  pitch,  roll,  and  yaw  were  10  deg,  15  deg,  and  30  deg,  the  actual 
(Euler)  yaw  angle  would  be  32.4  deg.  If  the  measured  pitch,  roll,  and  yaw  were 
30  deg,  40  deg,  and  30  deg,  the  actual  (Euler)  yaw  angle  would  be  47.0  deg. 

It  should  be  pointed  out  that  DTRC  has  been  fortunate  in  the  past  that 
most  measurements  made  on  the  carriages  and  at  sea  are  made  on  mild 
maneuvers.  In  the  case  of  radical  maneuvers,  such  as  the  emergency  recovery 
described,  the  section  of  the  run  that  is  of  the  most  interest  occurs  when  the 
pitch  and  roll  angles  are  small.  Therefore,  the  Euler  angle  assumption  is  a 


relatively  good  assumption.  However,  as  more  sophisticated  analytical  tech¬ 
niques  are  introduced,  the  error  in  yaw  can  be  important,  since  the  correlation 
of  body-centered  measurements  and  inertial  measurements  relies  on  the  coordinate 
transformation  matrix  which  uses  yaw  as  a  parameter.  For  example,  if  Eq.  4 
were  used  to  obtain  an  estimate  of  r,  the  used  in  this  equation  should  be 
obtained  by  converting  the  measure  yaw  to  Euler  yaw  by  using  Eq.  8.  Assuming 
that  the  measured  yaw  is  the  Euler  yaw  can  result  in  errors  in  r  on  the  order  of 
10*?^o  to  20%  for  radical  maneuvers. 

RADIO-CONTROLLED  MODEL  GYRO  PACKAGE 

The  radio-controlled  model  gyros  are  installed  in  the  model  in  such  a 
fashion  that  the  definitions  presented  in  Table  2  apply. 


Table  2. 

Position  vectors,  measurement  planes,  and  reference 
planes  for  the  radio-controlled  model  gyros. 

Measurement 

Position  Vector 

Measurement 

Plane 

Reference 

Plane 

Yaw 

Body’s  Y-axis 

Body’s  X-Y  plane 

Fixed  Y-Z  plane 

Pitch 

Body’s  X-axis 

Body's  X-Z  plane 

Fixed  X-Y  plane 

Roll 

Body’s  Y-axis 

Vertical  plane 

Fixed  X-Y  plane 

passing  through 

the  body’s  Y-axis 

A  comparison  between  the  Table  2  definitions  and  the  analogous  Euler 
definitions  (Table  1)  shows  that  the  roll  and  pitch  measurement  plane  definitions 
are  different.  This  difference  is  due  to  the  fact  that,  due  to  space  limitations, 
the  vertical  gyro  had  to  be  mounted  such  that  the  inner  gimbal  measured  roll 
and  the  outer  gimbal  measured  pitch.  To  use  this  gyro  measurement 
complement,  a  different  coordinate  rotation  sequence  must  be  performed. 

To  rotate  the  inertial  (or  fixed)  coordinate  system  into  the  moving  (or  body) 
coordinate  system,  the  first  rotation,  around  the  fixed  Z-axis,  rotates  the  fixed 
Y-axiz  under  the  body’s  Y-axis  (see  Fig.  3).  The  second  rotation,  around  the 
intermediate  X-axis,  rotates  the  intermediate  Y-axis  into  the  body’s  Y-axis  (see 
Fig.  3).  The  last  rotation  is  around  the  body’s  Y-axis  and  rotates  the  second 
intermediate  X-axis  into  the  body’s  X-axis  (see  Fig.  3).  Note  that  the  second 
rotation  is  measured  in  a  vertical  plane  and  the  last  rotation  is  measured  in  the 
body  X-Z  plane.  These  rotations  are  in  the  measurement  planes  of  the  gyro; 
therefore,  the  angles  measured  by  the  gyro  package  for  pitch  and  roll  will 
correspond  to  the  transformation  angles.  The  yaw  rotation  (the  first  rotation  is 
in  the  fixed  X-Y  plane  and  is  not  in  the  measurement  plane  of  the  yaw  gyro. 

The  body-to-fixed  coordinate  transformation,  described  by  the  above  rotations, 
is  shown  in  Eq.  9. 
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ip  (cosv;,cos0,  -  sin»^,sinf  ,sin0, 

jp  =  (sinv;,cos0, +  cosi^>,sini^,sin0,) 
kp  -cos<^,sin0, 


—  sim^jCosi^i 
COSJp,COSiti, 

sinif, 


(cosvj,sin0i  +  sinifi,sinf  ,cos0,) 
(sinif»,sin0,  -cosv;|Sin^,cos0,) 
COS<|i|COS0, 


'b 

Jb  (9) 


Where  ip,  jp,  and  kp  are  the  unit  vectors  in  the  fixed  system,  and  ig,  jg,  and  kg 
are  the  unit  vectors  in  the  body’s  system. 

The  variables  in  this  transformation  are  rotation  angles  for  the  new 
coordinate  transformation.  Because  these  angles  are  different  from  the  Euler 
transformation  angles,  the  subscript  I  is  used  to  differentiate  the  measurements 
from  the  Euler  measurements.  The  coordinate  transformation  should  produce 
consistent  results  when  translating  a  vector  from  one  coordinate  system  (fixed  or 
moving)  to  another  coordinate  system  (moving  or  fixed)  no  matter  what 
coordinate  transformation  is  used.  The  transformation  equations  to  transform 
the  RCM  measured  pitch  and  roll  can  be  obtained  by  equating  the 
corresponding  [3,1]  and  [3,2]  elements  of  the  Euler  and  the  RCM 
transformation  matrices.  These  transformation  equations  for  RCM  pitch  and 
roll  (RCM  to  Euler)  are  given  in  Eqs.  10  and  1 1 . 


6  =  sin  *(cos<t>,sinej)  (10) 

<ti  =  sin” *(sin4'j/cos6)  Note  0  from  above  equation.  (11) 

Where  and  <|)j  are  the  measured  pitch  and  roll,  and  0  and  are  the  Euler 
pitch  and  roll  angles. 

Yaw  presents  a  special  problem  since  the  Euler  and  the  RCM 
transformations  (Eqs.  10  and  11)  both  assume  that  the  measurement  plane  is  the 
fixed  X-Y  plane.  The  actual  RCM  yaw  measurement  and  the  Euler  angle 
measurement  transformation  equation  were  derived  by  transforming  two  vectors 
in  the  fixed  X-Y  plane  to  form  an  Euler  angle  yj  to  the  body’s  coordinate 
system.  The  angle  formed  by  the  projections  of  the  two  vectors  in  the  body’s 
X-Y  plane  is  the  angle  t^2  ^^at  would  be  measured  by  the  RCM  instrumentation. 
This  expression  was  simplified  to  obtain  the  transformation  equation. 

As  mentioned  earlier,  the  spin  axis  of  the  yaw  gyro  is  aligned  with  the  fixed 
X-axis.  The  measurement  plane  is  the  body’s  X-Y  plane  because  the  outer 
gimbal  is  being  used  to  measure  the  yaw  angle.  A  vector  in  the  intersection  of 
the  reference  and  measurement  planes  can  be  found  by  forming  the  vector  cross 
product  between  a  normal  to  the  reference  plane  and  a  normal  to  the 
measurement  plane.  The  angle  measured  by  the  gyro  is  the  angle  between  this 
intersection  vector  and  the  position  vector.  For  the  RCM,  the  normal  to  the 
reference  plane  is  the  fixed  X-axis;  and  the  normal  to  the  measurement  plane  is 
the  body  Z-axis.  Therefore,  a  vector  in  the  intersection  between  the  reference 
and  measurement  planes  is  given  by  Eq.  12. 

V  =  (T*ip  X  kg)  (12) 

Where  V  is  the  vector  in  the  intersection  between  the  reference  and  measurement 
planes,  and  T  is  the  Euler  fixed-to-body  transformation  matrix  (the  inverse  or 
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transpose  of  the  matrix  given  in  Eq.  1).  The  unit  vector  ip  is  the  unit  vector  for 
the  X-axis  in  the  fixed  coordinate  system  and  kg  is  the  unit  vector  for  the  Z-axis 
in  the  body  coorainate  system. 

The  Tangent  of  the  angle  between  the  vector  V  and  the  position  vector  is 
given  by  Eq.  13. 

tan  yj2  =  (|V  x  jD/fV-jg)  (13) 

Where  V  is  the  vector  calculated  above,  and  jg  is  the  unit  vector  for  the  body’s 
Y-axis.  Note,  this  is  the  magnitude  of  the  vector  cross  product  divided  by  the 
vector  dot  product. 

Solving  Eq.  13  for  the  Euler  yaw  angle  we  obtain  the  final  transformation 
Eq.  14. 

tp  =  tan“'  ((tantp2COsfl -t- sin<t>sin6)/cosi^).  (14) 

The  variables  0,  4.  and  in  this  equation  are  Euler  pitch,  roll,  and  the  yaw, 
where  xpj  is  the  actual  measured  yaw  measurement.  It  should  be  noted  that  Eq. 
14  provides  a  method  of  obtaining  the  Euler  yaw  angle  from  a  free  gyro  whose 
spin  axis  is  initially  aligned  with  the  body’s  X-axis.  This  transformation  equation 
will  be  different  if  the  spin  axis  is  aligned  along  the  model’s  Y-axis. 

The  error  associated  with  the  assumption  that  the  RCM  gyros  measure 
Euler  angles  can  be  seen  from  the  following  examples.  If  the  RCM  gyro  package 
measured  30  deg,  40  deg,  and  30  deg  for  pitch,  roll,  and  yaw,  respectively,  then 
the  Euler  angles  would  be  22.5  deg,  44.1  deg,  and  48.1  deg.  This  example  is 
typical  of  the  measurements  obtained  during  a  radical  emergency  recovery.  A 
comparison  of  the  yaw  angle  obtained  using  Eq.  8  for  the  same  gyro  output 
measurements  indicates  a  slightly  different  yaw  angle.  This  is  due  to  the  fact 
that  the  pitch  and  roll  angles  measured  by  the  RCM  gyros  are  not  Euler  angles 
due  to  the  necessity  of  mounting  the  pitch-roll  gyro  with  pitch  being  the  outer 
gimbal. 

If  the  RCM  gyro  package  measured  10  deg,  15  deg,  and  30  deg  for  pitch, 
roll,  and  yaw,  then  the  Euler  angles  would  be  9.7  deg,  15.2  deg,  and  32.4  deg. 
This  example  is  typical  of  the  gyro  data  obtained  from  a  mild  emergency 
recovery. 

We  should  point  out  that  the  RCM  gyros  have  accuracies  on  the  order  of 
0.25  deg.  This  implies  that,  for  ail  mild  maneuvers,  the  measurement  error 
associated  with  the  gyros  is  larger  than  the  errors  that  are  experienced  by 
assuming  that  the  gyros  measure  the  Euler  angles.  It  should  be  pointed  out  that 
a  great  majority  of  the  maneuvers  performed  in  the  past  have  been  “mild 
maneuvers.’’ 

The  RCM  data  are  presently  processed  by  a  Kalman  filter.  We  noted  that 
the  Kalman  filter  results  are  consistent  (white  residuals)  for  the  mild  maneuvers. 
However,  the  Kalman  filter  residuals  for  radical  maneuvers  are  not  white.  The 
non-white  residuals  have  required  a  case  by  case  analysis  procedure  to  derive 
time  histories  of  the  state  variables  to  use  as  inputs  for  the  System  Identification 
process.  This  problem  is  due,  at  least  in  part,  to  the  assumption  that  the  RCM 
gyros  measure  Euler  angles. 
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CONCLUSIONS 


The  RCM  gyro  measurements  for  pitch,  roll,  and  yaw  can  be  converted  to 
Euler  angles  by  using  the  transformation  equations,  9  and  13,  that  are  reiterated 
here  for  convenience. 


Pitch 

6  =  sin“*  (cos<t»jsin0i). 

(10) 

Roll 

iji  =  sin~*  (simfij/cosS),  (note  9  from  Eq.  10) 

(11) 

Yaw 

ip  -  tan  “  ^  (tantp2cos0  -i-  sin<^sin0)/cos(|>)  (note  6  and  from 

Eqs.  10  and  11). 

(14) 

Here  t^2>  measured  yaw,  pitch,  and  roll,  and  6,  and  are 

the  Euler  yaw,  pitch,  and  roll. 

Note  that  these  equations  can  be  used  to  convert  gyro  measurement  angles 
to  Euler  angles  as  long  as  the  gyro  instrumentation  consists  of  a  vertical  gyro 
with  the  outer  gimbal  mounted  in  pitch  and  a  free  gyro  (for  yaw)  with  the  outer 
gimbal  mounted  to  measure  yaw.  Because  this  gyro  complement  is  relatively 
inexpensive  and  readily  available  when  compared  to  the  cost  of  a  stable-table 
system,  these  transformation  equations  can  result  in  considerable  cost  and  time 
savings  for  model  tests. 

In  conclusion,  the  results  of  all  DTRC  experiments  would  be  improved  if 
the  transformation  equations  presented  in  this  report  were  applied  to  the  gyro 
measurement  data. 
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APPENDIX  A 

SUPPORTING  CALCULATIONS 

To  transform  vector  quantities  from  the  model  coordinate  system  to  the 
fixed  coordinate  system,  the  below  rotations  are  required. 

1.  Rotate  around  the  model  Y-axis  until  the  model  X-axis  is  in  the  fixed 
X-Y  plane. 

2.  Rotate  around  the  X-axis  until  the  Y-axis  is  in  the  fixed  X-Y  plane.  (The 
X  and  Y  axes  referenced  here  refer  to  the  model  coordinate  system  after 
the  rotation  described  in  step  1  has  been  performed.) 

3.  Rotate  around  the  Z-axis  until  the  Y-Z  plane  aligns  with  the  fixed  Y-Z 
plane.  (The  Z-axis  and  the  Y-Z  plane,  reference  here,  refer  to  the  model 
coordinate  system  after  the  rotations  in  steps  1  and  2  have  been  performed.) 

These  rotations  can  be  described  in  matrix  notation  as  Eqs.  A.l  through  A. 3. 


cosS, 

0 

sindj 


0 

1 

0 


-  sin0j 
0 

cos8j 


'b 

Jb 

*^B 


‘T 

i" 

h 

k" 


1 

0 

0 


0 

COS<frj 

-sin^, 


0 

sin4, 

COSI^, 


kj 


‘f 

h 


costal 

-sini/^j 

0 


sini^j 

COStf^j 

0 


0 

0 

1 


Jt 


(A.l) 


(A.2) 


(A.3) 


Here  the  rotation  angles  are  0j,  <^j,  and  The  unit  vectors  before  the  rotations  are 
“B’  ~T*  “T  right-hand  side).  The  unit  vectors  after  the  rotations 

are  -j,  -j,  and  -p  (on  the  left-hand  side).  An  examination  of  the  columns  of 
any  of  these  rotation  matrices  shows  that  all  columns  have  a  magnitude  (square 
root  of  the  sum  of  the  squares  of  the  elements)  of  1 ;  and  that  the  dot  product 
of  any  column  with  either  of  the  two  remaining  columns  is  0.  This  implies  that 
the  inverse  of  each  of  the  individual  matrices  can  be  found  by  taking  its 
transpose.  Because  any  transformation  matrix  is  composed  of  the  product  of  the 
individual  rotation  matrices,  the  inverse  of  the  transformation  matrix  can  be 
found  by  taking  its  transpose. 

The  transformation  matrix  to  transform  vector  quantities  from  the  model 
coordinate  system  to  the  fixed  coordinate  system  can  be  obtained  by  multiplying 
the  above  rotation  matricies  together  in  the  order  that  the  rotations  were 
performed.  This  multiplication  yields  Eq.  A. 4. 
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ip  (cosv/|COs6| -sini^,sin^,sin0,)  -sinjf;,cosi^,  (cosy, sinflj +sin:^,sin^, cose,)  13 

jp  =  (cosi<>,sini^,sine, +sinviCOse,)  costp,cosi^,  (sinv,sin0, -cos4;,sin^|Cose,)  j0  (A.4) 

kp  -cosi^isine,  sin^,  cos+[Cos0,  kg 

Where  0p  itij,  and  ^p^  are  the  rotation  angles,  ip,  jp,  and  kp  are  the  fixed  unit 
vectors,  and  ig,  jg,  and  kg  are  ne  unit  vectors  in  the  model  coordinate  system. 

The  rotation  angles  0,  and  <}>,  are  tuc  ngles  measured  by  the  RCM  gyros. 
For  yaw,  the  RCM  gyros  measure  the  angle  in  the  body’s  X-Y  plane  rather  than 
the  required  fixed  X-Y  plane.  It  is,  therefore,  necessary  to  translate  the 
measured  yaw  angle,  t^2>  desired  yaw  rotation  angle,  This 

transformation  will  be  derived  by  calculating  the  angle  between  the  intersection 
of  the  measurement  plane  and  the  reference  plane;  and  the  position  vector  (the 
model’s  Y-axis). 

For  the  yaw  gyro,  the  spin  axis  is  aligned  with  the  fixed  X-axis.  Therefore, 
the  fixed  X-axis  is  the  normal  vector  to  the  yaw  reference  plane.  Because  the 
outer  gimbal  is  being  used  to  measure  the  yaw  angle,  the  model’s  X-Y  plane  is 
the  measurement  plane.  A  vector,  V,  in  the  intersection  of  the  reference  plane 
and  the  measurement  plane  can  be  obtained  by  taking  the  vector  cross  product 
of  the  normal  vectors  to  the  planes  (Eq.  A. 5). 

V  =  Tip  X  kg  '  (A.5) 

Here,  V  is  the  vector  that  is  in  the  intersection,  and  T  is  the  transformation 
matrix  from  the  fixed  coordinate  system  to  the  model  coordinate  system  (the 
transpose  of  Eq.  A.4).  The  vector  ip  is  the  unit  vector  for  the  X-axis  in  the 
fixed  coordinate  system  and  kg  is  the  unit  vector  for  the  Z-axis  in  the  model’s 
coordinate  system.  Using  the  definition  of  the  cross  product,  V  can  be  found  by 
taking  the  determinant  of  Eq.  A. 6. 

ifl  is 

(costf^icos0, -sini^isiniji,sin0,)  -sint^, cos'll  (-f-costf»,sin0, +sinif(,sin4,cos0,)  (A.6) 

0  0  1 

Computing  this  determinant  yields  Eq.  A.7. 

V  =  -sinv^jCOs<t>jig  +  (-cosi^^jCOsSj +  sintf»jsin<tijsinej)  jg  (A. 7) 

The  RCM  yaw  gyro  measures  the  angle  between  the  vector  V  and  the  model’s 
Y-axis,  jg.  The  tangent  of  this  angle  can  be  found  by  dividing  the  magnitude  of 
the  vector  cross  product  of  V  and  jg  by  the  vector  dot  product  of  V  and  jg. 
Performing  this  operation  gives  Eq.  A.8. 

tanj^2  =  cos<|>j/(cotif;jcos0, -sin^jSinSj)  (A. 8) 

The  expression  relates  the  measured  yaw  angle  to  the  rotation  angle  given  in  Eq. 
A. 3.  This  derivation  is  included  because  noise  correlation  considerations  may 
dictate  that  it  would  be  better  to  use  the  gyro  measurements  for  pitch  and  roll 
rather  than  to  transform  the  measurements  to  Euler  angles. 
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The  measured  angle,  ip2>  related  to  the  Euler  angle  by  using  the 

Euler  transformation  in  the  computation  of  the  vector  V.  When  this 
computation  is  performed,  the  vector  V  is  given  in  Eq.  A. 9. 

V  =  ( -  sint^cos*}*  +  sin<t»cosipsine)i3  -  cosipcosSjg  (A. 9) 

The  expression  transforming  the  RCM  gyro  yaw  measurement  to  the  Euler  yaw 
angle  can  be  found  by  dividing  the  magnitude  of  the  vector  cross  product  of  V 
and  the  model  jg  unit  vector  by  the  dot  product  of  V  and  the  model  jg  unit 
vector.  The  resulting  expression  is  then  solved  for  and  results  in  the 
transformation  equation  given  in  Eq.  A.  10. 

yj  =  tan“*  [(tantp2cos0  +  sin<|>sin6)/cosi^]  (A.  10) 

Where  9  is  Euler  pitch,  <j)  is  Euler  roll,  ip  is  Euler  yaw,  and  1^2  is  the  measured 
yaw  angle  from  the  RCM  gyro  package. 

The  equations  that  transform  the  RCM  pitch  and  roll  measurements  to  the 
Euler  angles  can  be  determined  by  equating  the  Euler  transformation  (body  to 
fixed)  with  the  transformation  given  in  Eq.  A. 9.  The  transformation  matrix 
from  the  body’s  coordinate  system  to  the  fixed  coordinate  system  to  the  fixed 
coordinate  system  should  be  unique  no  matter  how  the  rotation  angles  are 
measured.  Working  with  the  [3,1]  elements  for  pitch  and  roll,  we  obtain  the 
translation  equations: 

6  =  sin"*  (cos^jsinflj)  (A.  11) 

=  sin"*  (sin<^j/cos0)  (A.  12) 

using  9  from  Eq.  A.ll. 

This  technique  could  also  be  used  to  relate  ipj  to  the  Euler  yaw  angle; 
however,  tpj  is  not  the  angle  measured  by  the  RCM  gyros  because  it  is  measured 
in  the  fixed  X-Y  plane  and  the  RCM  gyros  measures  the  yaw  angle  in  the 
model’s  X-Y  plane.  This  would  require  another  transformation  equation. 

Because  we  have  no  way  of  measuring  the  translation  equation  was  not 
derived. 

As  mentioned  earlier,  the  yaw  angle  being  used  in  the  Euler  transformation 
is  measured  in  the  fixed  X-Y  plane.  This  implies  that  the  direct  measurement  of 
this  angle  will  require  a  stable-table  to  keep  the  measurement  plane  of  yaw  in 
the  fixed  X-Y  plane. 
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